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A RANDOM DISTRIBUTION OF RADIAL MOTIONS 


William R. Thicks tun, Jr. 

Theoretical Mechanics Branch 

ABSTRACT 

Two quite unrelated problems motivated the efforts described in this 
report. One pertains to the question of whether a singularity in space- 
time can develop from an initially non-singular distribution of matter. 
The other problem pertains to the initial distribution of comet orbits 
about the sun from which certain "capture" theories attempt to predict 
the observed number of short period comets . For both of these problems 
it is instructive to start with a model consisting of a random distrubu- 
tion of particles constrained, however, to strictly radial motion in a 
central gravitational field. 


iii 


1 . 


INTRODUCTION 


Imagine a swarm of particles with random velocities and distances 
from a suitable origin. For example a swarm of comets about the sun, or 
a cluster of galaxies under their mutual gravitational attraction. 

Further suppose the system is in quasi-equilibriun, i .e . only very slowly 
permit the density of the distribution to change, e. g. because of random 
close encounters or a weak resist ing medium. It is instructive to con- 
sider a grossly simplified version of such a swarm of particles by 
restricting attention to purely radial motions. This model may be moti- 
vated by two quite unrelated problems which we may take up in turn. 

First, for a swarm of comets it should be noticed that we on Earth 
OfJ y oee or observe those that come within one or two Astronomical 
Units (A.U.) of the Sun. Hence that sub-set of the set of all long 
period comets contains only those whose orbits are long slender ellipses 
with eccentricities close to one. Thus we only observe those that from 
time to time drop almost radially toward the Sun from a great distance. 
One might at first glance get the impression that such comets are sort 
of kept out in "deep storage" for centuries at a time only to drop in 
for short visits. In a sense this is true as we shall show. Further, 
by varying the randomness of the distribution one might hope to throw 
some light on an old question about the short period comets usually 
attributed to Laplace but most recently studied by Everhard^ . 

Secondly the model to be discussed in this report may be reviewed 
as an attempt to start with a simple easily visualized Newtonian 
picture and then by slowly increasing the density at the center try, 
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so to speak, to creep up upon a gravitational collapse as predicted by 

general relativity. In the first section of this report then we shall 

discuss the random distribution of radial motions in Newtonian Mechanics. 

Then we shall review briefly the classical solution of Oppenheimer and 

(2) 

Snyder on gravitational collapse. Finally we shall view collapse 
from the point of view of a non-uniform density distribution. 

2 . THE MOTIONS IN NEWTONIAN MECHANICS ABOUT A CENTRAL BODY 

Imagine the swarm of particles with random velocities and distances 
from the origin. We desire the?'r distribution for all subsequent time. 
The equation of motion for each is given by 

ci^x* M 

(2.1) = - G — , (M = total mass, G the grav. const., r the 

dt 2 r 2 distance) 

if one assumes each can be regarded as moving in a central fj M and 
ignore the shielding effect which occurs when the motion cf the one 
particle is within most of the mass of the swarm. 

The motion of each particle is next regarded as purely radial. This 
idealizes the motion of a long slender ellipse. The motion may be des- 
cribed as follows. (These may be verified by differentiating, or by 

( 5 ) \ 

referring to a standard text such as Szebehely'*' .) 

(2.2) C 3 (1 + cos f) 

t' = U + sin *), 
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where the new radial variable £ instead of r is 


(2.3) C - r/a 

and a new time t ' instead of t is 

(2.1,) f . jsL t 

In these variables the nation is in the sense of increasing Equation 
(2.2)furnishes, for i(r = 0, £ = 2 , and t' = 0. Subsequently £ decreases 
to zero as i|r increases to tt; for which value t ' = tt . The plot of 
the motion in the (£, t') -plane is a cycloid, with its cusp cor- 
responding to infinite speed at the origin ( occurring at t 1 = tt ) . 

The particle then rises back to its maximum distance at £ = 2 , or in 
the original variable, at r 3 2a * Denote this value by t. 

Before leaving these familiar relations, observe that the single 
particle spends most of its time out near r * t and very little near the 
origin. It will be worth stating this precisely in what follows by 
defining a density distribution for a single particle whose integral 
over the whole radial path is the mass m of the particle. 

To do this observe either by going back to (2.1) and obtaining the 
first integral or more directly from (2.2) that 



(2.5) 



hence the reciprocal of 


dt 


, is proportional to the relative time spent in 
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the interval Q to Q + d£ . Let us denote the density by cr^, with 


( 2 . 6 ) 


P 

j o^(r) dr = m > the particle mass. 

o 


Then one has 


(2.7) a/r).|-a 

1 }n 1 „K3jjn 


1/2 


— i 


v" 

vr* 


% 




Notice , as stated, that ct. (o) = 0, whereas a (l) is infinite, so indeed 

V { 

in this sense such a particle is indeed out in "cold storage" and only 
"drops in" from time to time. 

If now we have a swarm of particles n^ in number, each of imss m, 
all in radial orbits of major axis £, but uniformly distributed in all 
directions, we may define a density p (r) for all of them defined so that 

*V 

the integral over tie whole space out to r = .(, will be n m: 


( 2 . 8 ) 


n.m 


,(r) Aik 


( r A) 


1/2 




, 0 < r < l 


, r i» l 


Notice p (o) is no longer zero because of the large supply from 

v 

all directions that enter a small neighborhood of the origin. 

Let us next make the swarm less artificial by permitting a wide 
distribution of energies, characterized by t, to be present. Indeed 
let us introduce a distribution in l, namely jx(4) with now being the 
number with maximum distances from t to l + dt . Then for the density 
distribution p(r) for the whole cluster, we have 
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where the lower limit moat be r, not zero since the function p t (r) 
vanishes for t £ r. Thus 


(2.10) p(r) = -2- r' 3/2 f . 

6tt 2 r t ^-r 

Before continuing we should discuss the consistency of this equation. 
It appears to define p(r) for a given distribution of orbit sizes ^(-t) . 
However we assumed the central force law (2.1) in the derivation. Ufcus 

(2.10) is correct for the density of comets, for example, whose total mass 
is negligible compared to a central body such as the sun; but (2.10) 
will not hold for a swarm of particles in their mutual field because of 
the shielding effect. We shall discuss this latter point in the next 
section . 

To get some feeling for (2.10) and to keep the integrations easy, 
let us consider a simple power distribution for ^(-t)* and write 

(2.11) n(t) = cl, 1 " 11 . 

If we let N denote the total number of all energies out to some 
remote last value of l say L, we may chose c so that 

L 

J o iiU)dt - H. 


( 2 . 12 ) 


r l u U- r) A/ * 


° (l+x)“x x/t 


or as we prefer to write 


(2.14) 


r( n - ■§) i _ 

^•Tra— ^ r * 


Thus from (2.10) we nave 


( 2 .1 5 ) p(r) = C n r- (nn) , 


Finally integrating over a sphere of radius r we find for the mass 


(2.16) H(r) • C r 2 ' n . 

.;-n n 


To sum up as follows: For any arbitrary distribution n({,) of 

orbit sizes the corresponding density distribution p(r) is given by (2.10) 


In the special case of all having the same size t the density is given 
by (2.8). In the special case of ji(-t) given by a plausible power law (2.11) 
the density is given by (gjjj) and the mass interior to a sphere of radius 
r is given by ( 2 J6) . Thus n=»2 corresponds to the usual infinite raass-at- 
the-origin, central field approximation . Cases with n < 2 correspond to 
less concentration, and indeed n = -1 gives a uniform density. 

Example : Suppose an astronomer is aware at any one time of 5 comets in 

our neighborhood (i.e. a cube 1 A.IJ. on a side centered at r - l) whose 
major axes are known to be 1000 A.U. How many are there altogether? 

Ans. 10,000,000. 


3. THE MUTUAL MOTIONS IN NEWTONIAN MECHANICS. 

As mentioned before, to treat the case where there is no large 
central mass bo that particles in the inner regions of the swarm are 
shielded, we must give up equation (2.1). To correctly account for the 
shielding even in this setting of pure radial motion is difficult, it 
involves integrating r » - V MfrJ/r 2 for M(r) arbitrary, working through 
to an equation analogous to (216) and then regarding tne latter as an 
inte ,ral equation for that n(<,) which furnishes the given M(r) . We 
shall not attempt this in general but consider the problem for M(r) 
given as a power law in the same spirit as the last section. 

Indeed for 


(3-1) 


c -( 1+n) 

C n r 
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we have 


P ** 

(3.2) M(r) = j pr 2 dr dfi =» 4rr J pr 2 dr 
V o 


, r _ 1-n , 4tt C n 2-n ^ _ 

= 4tt 1 C_ r dr =» r , n < 2 


J n 

o 


2-n 


and hence 


u 4 tt C G 
•• „ M n -n 

r = - G — =» - — r 

2 2-n 

r 


Let 


4tt C G 

* = 2^r > 30 


(3.3) 


r = -ft r 


-n 


§ 


- J E - r 1 " 11 , n / 1 


J E - 2 -ft In r , n = 1. 


VThere E is iasrely an integration constant . 

Finally t may be obtained by quadrature as an incomplete gamma 
function. However the behavior of the solution is easily found from 
(3-4). The principle effect of the shielding is to remove the infinite 
velocity at the origin for n < 1. If ve restrict attention to bounded 
orbits and choose E so that r(o) * l, the motions qualitatively may be 
sketched as in the Figure 1. 





FREE FALL PATHS 


Figure 1 


Notice that n = 1, the logarithmic case, serves as the dividing line 
between those orbits that return on themselves and those that cross the 
origin (The understanding here is tnat the radial paths are thought of 
as approximating long slender ellipses in the n= 2 case) . The n ^ - 1 
case is that of a uniform dust cloud of constant density. 

To continue as described after equation (2.5 ) > we have from (3.4) 


(3-5) ^(r) = 




where a is a suitable constant chosen so that the integral of a 
taken over r is the mass of the particle . We eschew the detail here of 
explicitly expressing a in terms of m and n by suitable gamma functions 


If wf: pursue the more general case through to the analogue of equation 
(?. .10) we obtain 


(5.6) p(r ) , « 4:° ..-. f 4iL^_ 

4nr 2 r l-n^l-n 


or introducing x through l = rx 


(3-7) p(r) 


Of aA-p r s 

4tt 


r u (rx) dx n- 

/— ' S ~~2 

1 V x 1 " - 1 


Now this relation is based on the assumption (3 *1 ) , hence we seek a 
function for which this is the case. But from the form of ( 3 . 7 ) 

it is clear that hK) must be a power function in l if p is to be a 
nower function in r. Thus we assume 


(3-8) 


M. = , 


where m is an as yet unknown power. Absorbing all the constants and 
equating ( 3 .l) to (3*7); 


C r- (1+n) = o' „ 


hence we have a solution provided 


(3-9) 


C = C , 
n n* 



(4.0) 


1 3 
m = 2 ' 2 


Thus we have a family of solutions to (5*7) depending upon n. 

The functional dependence of the va-lous variables may be summarized in 
a table as shown in Table 1. The extreme central concentration n = 2 is 
on the left, the complete lack of concentration is n * -1 on the right. 


variable 

power 
of r 

n»2 

n=»l 


n=0 

n= -1 

P 

-1-n 

-3 

-2 

1-33 

-1 

0 

r 

-n 

-2 

-1 

-0.33 

0 

1 

M- 

2 2 

-2.5 

-1 

0 

0.5 

2 

E-r 2 

1-n 

-1 

( logr) 

0.66 

1 

2 

M 

2-n 

0 

1 

1.66 

2 

3 


TABLE 1 


The case n = - has a rather nice interpretation. Consider a system of 
particles initially at rest and uniformly distributed and suddenly 
allowed to fall under their mutual gravitation. Tien if we allow their 
energies to be unchanged, but the phases of their motions to drift 
incommensurably, the state of motion can shift statistically from n = -1 
to n = l/3- It may be of interest to pursue the question implied here 


more seriously in view of the recent efforts of Barnes and Whit row 


(19) 


and Penston 


( 20 ) 


However we shall not follow this line of argument 


further in this report 



h . RELATIVISTIC GRAVITATIONAL CONTRACTION 

In their now classic report on gravitational contraction (1929), 

(2) 

Oppenheiraer and Snyder obtained a solution of the Einstein equations 
of general relativity that correspond to the fall of a ball of matter of 
uniform density. Indeed it had previously been shown for sufficiently 
large mass concentrations that "cold" matter must contract^ . The 
implications of this, especially for stellar theory has prompted numerous 
papers since ^ . The literature is in fact much too extensive to 

review here. The basic conclusion seems inescapable and has led 
Geroeh^ , Hawking' 1 ^ , and Penrose^ 12 ^ to propose rigorous mathematical 
proofs, under fairly broad assumptions, that such collapse must occur. 

One feature of the original Oppenheimer-Snyder paper that is 

bothersome, is the fact that to obtain a first result of this kind they 

quite naturally considered the simplest case. It turns out that a 

uniform cloud of particles, Co-moving, will contract steadily and 

uniformly down to a point singularity. In spite of general theorems 

such as that of Penrose it is interesting to look at other specific 

( 13 ) 

examples. One recent example due to Israel considers a thin 
spherical shell of matter that collapsed due to its own gravity. 

Israel would have it rebound; appearing to go backwards in time after 
the "bounce". Other examples suggest themselves. Before going into 
them, we shall review the principle results of Oppenheiraer and Snyder 
here. They assume, for an observer at rest at a great distance, the 
Schwarzscliild metric given in terms of coordinates (r, 0, cp, t) by 


( 4 . 1 ) ds 2 =» (l - ^) dt 2 - (l - ^r )" 1 dr 2 - r 2 (d 0 2 + sin 2 © dtp 2 ) 


For the falling matter they introduce co-moving coordinates 


(R, 9, cp, t) with metric 


( 4 . 2 ) ds 2 = dT^ - e 10 dR 2 - e® ( d9 2 + sin 2 9 dcp 2 ) , 


where aj and ou are functions to be determined. Since the matter is at 
rest in these coordinates and any internal stresses and pressures arc 
taken as negligible, all components of the stress tensor vanish except 
T 4 = p(R) . Substitution into Einstein's equations then furnish equa- 
tions for cu and u)i indeed they find 


(4-3) 


e 10 - e® cu tZ /ki?( R) , 


where f(R) is arbitrary. A restricted class of solutions is found for 
f(R) =1 for which they find 


(4.4) 


e* * (Ft + G) 4 ^ 3 , 


with F and G arbitrary functions of R. Since the equations are invariant 
under a coordinate transformation, such as taking R to be a function of a 
new variable R*, only one of F and G is really arbitrary so they set 
G = R 3 / 2 . If now the density p(R) be given initially as po(R) they find 


(4.5) 


FF' =■ 9" R 2 p 0 (R) . 
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We have repeated much of the argument here because we wish to use these 
equations for p 0 (R) chosen as in Section 3- \t this point the easiest 
case is to take p 0 (R) =* const. However such a distribution cannot ex- 
tend to all R but must fall off faster at infinity; so they introduce 
a value such that po(R) = const. R s, R^, p 0 (R) = 0 for R > R^. One 
then wishes to observe the collapse in the (r, 8, cp, t) coordinates. 

The transformation is obtained by finding r = t(R,t), t = t(R,T) using 
( 4 . 3 ) and (4.4) to transform (4.2) over to (4.1). Skipring the details 
here, we shall simply give the result ( r 0 being the Schwarzschild radius 
corresponding to the total mass.): 

R (1 - | ro^" 3 / 2 t)" 7 ° ; R < R fe 

R (1 - | r 0 ^ R -3 / 2 t ) ; \ 

f o / / • r 8 + r 0 2 

(R 3/2 - r 3 ' 2 ) - 2(rr 0 )* + r 0 In -j t, R > R. 

\r 0 2 r* - r 0 ® 

(4.7) t 4 

/M(y), R<R b 



i 


• J 
' 


v *"■ 1 


- : " • 


j 

V ’ ’ 


« 




where 

(4.8) 


y = \ [(f^) - 1] + Rj^r/roR. 


Here M(y) must be chosen to satisfy continuity at R = which furnisher 


(4.9) M(y) = | r 0 ”^( Rjj 3 ^ 2 - ro^y 3 ^) • 2r 0 y^ + r 0 tn . 


The appearance of the solution is easiest to visualize by appropriate 
graphs. Figure 2, gives the paths R = const, in the mixed (t,t) -plane. 
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Figure 3 shows the motion in (r,t). The latter shows the phenomena 
known as cut off. If we denote the moment of complete collapse by t 

CR 

we see that the falling motions appear to slow down for t in the 

neighborhood of t and rather quickly approach asymptotic distances 
UR 

which are the Schwarzschild radii for the total mass interior to each. 

That is tim =» r 0 , £im f ^ = T o> etc. each finding a limiting posi- 
tion. Since the asymptotic values are in fact approached rather quickly 
for t > t-, p one n&y ask if this is an effect of the assumption that all 

oil 

the particles started out together and hence all converged at the same 
time. We may also wonder at the consistency of this picture with the 
Newtonian picture of Section 3 where particles may fall from different 
initial points at different times. 

To answer some of these questions let us first point out an 
important distinction between the particles falling in Section 3 and 
the co-moving coordinate R following the motion here . Consider a small 
volume of space in the neighborhood Y, of a point at distance R initially. 
The volume f at any instant contains mass. The momentum components T. 
of the stress tensor will be zero even for matter in motion, so long as 
the net flux into or out of Y is zero. Thus R, the "coordinate-following- 
the-motion", only has this meaning in the sense that the center of mass 
of Y has a motion described by r a r(R;t) such as shown in Figures 2 or 
3 . Thus individual particle motions different from tne mean collapse 
motion can be understood to be encompassed by the relativistic theory. 

On the other hand the picture in Section 3 suggests that we look at other 
motions following equation (4.5) with po(R) chosen from Table 1. Per- 
haps we can sort of sneak up on the singularity without the coraaon 
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\ 


suddeness of t a t shown in Figure 2. We shall indeed take this up 
Cn 

in Section 5* 

Before leaving the subject of uniform collapse, it is worth one 

more fact of note. If each region characterized by the coordinate R is 

furnished with a clock, the deck will appear to stop at a certain 

proper time t before the critical value of t . Moreover this occurs 
L 

asymptotically as R appears to sink down to its ultimate lowest value, 
say r T . This defines a curve r_ = r T (T_) in the (r,T) -plane beyond 
which further collapse is unobservable to our distant rest observer. 
This may be found from (4.7) (corresponding to y*l) , with a little 
algebra, and is shown on Figure 2 as the dashed line. 

5- COLLAPSE FOR ALTERNATIVE DENSITY DISTRIBUTIONS 

In this section we shall adopt an assumed power law 

(5.1) p * CR U 

where u * - (l+n) in terms of the exponent in (3*l)> and introduce the 
co- moving coordinate R «* r(0) in place of r(t) . Following Oppenheimer 
and Snyder, as described in Section 4, the only departure comes after 
(4.5). We shall follow their argument for general u. Actually looking 
at Figure 2, we see that by taking po(R) " const, for R < R^, po(R) * 0 
for R > they have in a sense covered two cases: u ■ 0 and u ■ -3. 

So that for 0 < u < -3 we may expect the paths of constant R in the 
figure to resemble those labeled v> ■ *3 rather than all converging to 
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j as in the u = 0 case. That this is indeed the more general situation 

CK 

we shall see shortly. 

Let us integrate (4.5) then to find 


18 tt Cu 

3 + v 


R 5 ^ ♦ 


constant. 


Dropping the non-essential constant we may write 


l5n Cu 

U 


an 

(5.2) F = -Y a 2 , Y 

hence e^ 2 = r furnishes in place of (4.6) 


/ Ion 

V 3 + 


(5-3) 


r = R(l- YtR L/2 ) 


Thus this motion is a solution of Einstein's equations and indeed should 

describe a possible relativistic behavior of the swarm of particles 

described in Section 2. Thus even though the motion is a "non-zero- 

temperature-motion", collapse can still occur. (This is in large part 

due to the zero angular momentum of pure radial motions.) To 6ee what 

the motion looks like to our distant observer we need, in addition to the 

function r = r(R, r) given by (5-3), a relation t ■ t(R,r). Tnis is 

( 2 ) . 

obtainable (still following ) by solving 

|t(R,r) |r(R.r) b &t/SR 

for t(R,r) regarded as a first order partial differential equation. 


IS 


since r*r/^T and *r/*R art obtained from (5.3)- Thus one finds 


(5-4) 


2 1+ T 

% V R 2 


(1-YtR U//2 ) 


V5 



(l-Y T R U//2 ) ‘ 5 

(1 -Y & R u/2 ) 


at 

an 


0. 


If <P 3 <p( R, t) be any first integral of this equation the general solution 
is t » \jr(ct)) > for if an arbitrary function. Note however that for all 
u < 0 the density given by (5-l) decreases to zero as R -• ®, so we wish 
the mapping (R,t) - (r,t) to approach the Schwarzschild field asymptoti- 
cally for all t . This fixes the function ^ < At this point one could, 
for given u, carry out the integration of (5*^) by setting up the 
equations 


(5.5) 


dR _ ^ 
dt 3 y 



L\ + ,.Y T . 


R U / g l B/3 



Y T R U//2 


If M 3 poCRjR 4 * dR be the total mass of the system, and we denote 

r the Schwarzschild radius corresponding to M, we may row write the 
s 

transformation as 


by 


R( 1 - Y t R u/2 ) 


2/3 


(5.6) 




2/3 l /2 

<p ♦ 2<p - tn 


JsuLl 

- 1 


where <p <p(R, T ) is any solution of (5.5) . 

Of particular interest will be the solution * 1 which 

describes the locus in the (R,t) space of all points that map to t * • 
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in the (r,t) space. These are the ultircate-lowest-apparent-radii 
reached. It is of passing interest to note that (5*5) is separable 
for two cases, y = 0 and u =* -2. Even in the u = -2 case, however, a 
numerical integration is necessary. 

Rather than get into a detailed solution of (5-5) we will content 
ourselves with a few comments on the nature of solutions . 

Figure 4 furnishes a sketch of the curves rar(r) described by (5-3) 
for v = 1> u = -1. These have been extended across the r=0 axis in the 
spirit of a "rebound" from the singularity. It should be noticed that 
two difficulties occur. For one thing if we compare with the Newtonian 
figure 1, we see that for all u, not only n s 0, one has dr/dj infinite. 

r 



FIGURE 4 
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This is a reflection of the fact that even though r is only a coor- 
dinate variable, none-the-less r=0 is the apparent location of a true 
singularity since the intrinsic curvatures are infinite . 

Secondly, as continued across r=*0, neighboring curves cross forming 
an envelope so that the coordinate transformation (R,t) -* (r,t) is no 
longer single-valued. This leads to two questions. One: what is the 

equation of the envelope? Two: Does the mapping ever break down with 

such a fold before the axis is crossed? The answer to the latter ques- 

u/2 u/2 

tion is no, since for u < 0, R 2 > Rj. implies Ra < R 1 and hence 


u/2 u/2 

1 - Y t R 2 >1-YtR , 

provided 1 - y t R 1 u ^ > 0. Thus rz(j) > ^(t) if ra(o) > r x (o) for all 
t for which r x > 0. 

That an envelope does in general exist for negative r is found from 
(5*3) by direct construction, i.e. eliminating R from (5-3) using 
dr/'cR = 0. The equation of the envelope is 


" <1 ♦ $>* ( ! * i? 


2/ 3 2 

($> T- U 


The case y * 0 degenerates to a point at t = t^ r , the case u = -3 gives 
r * V 2/3 t 2 ^ 3 . For u = -1 we have r proportional to t 2 as shown in 
Figure 4. 

The only really basic difference then between the case u * 0 and 
other power laws is that instead of the whole space seeming to simultaneously 
settle down to the ultimate state as in Figure 3, one has the inner regions 
effectively cut off first, followed by the outer regions later, so that 
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one must view a steadily expanding "black hole". 
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